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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2962. | 


ANALYTICAL INVESTIGATION ОЕ rcm LIMIT FOR | DIAMOND-SHAPED E 
ATRFOIL IN TRANSONIC AND SUPERSONIC "uw 


By Edmund Е. —— and. John 8. Serafini _ 


SUMMARY - 

Calculations have been made for the icing limit of a. diamond air- 
foil at zero angle of attack in terms of the stream Mach number, stream 
temperature, and pressure altitude. The icing limit is defined as a 
wetted-surface temperature of 229 F and is related to the stream condi- 
tions by the method of Hardy. | | 


The results show that the point most Likely to. ice on the airfoil 
lies immediately behind the зыш and is ирден ‘to possible icing 
at Mach numbers as high as 1.4. . 


їңористТОН a 


As the operational speed of aircraft .is increased through ‘the 
transonic region, the frictional heating available to prevent the forma- 
tion of ice on the aircraft becomes an important quantity. . A study of 
the probable icing limits of & high-speed airfoil in transonic and super- 
sonic flow was made in terms of the pertinent variables of flight Mach 
number and free-stream conditions of pressure and temperature. The pro- 
cedure used to determine the icing limit is based on the method presented 
in reference 1 for calculating the surface temperature of an insulated 
body (а body within which heat is not conducted from one section tó 
another) running fully wet in an air stream. The icing limit at any 
point on the airfoil was assumed to depend upon maintaining a tempera- 
ture of 320 Е at that point. Тһе free-stream static temperature corre- 
sponding to this icing limit was calculated for each point on the airfoil 
for particular values of free-stream Mach number ‘and pressure. 


The results presented herein were calculated for a symmetrical 
diamond airfoil at zero angle of attack for a range of airfoil- 
thickness ratios from 0.02 to 0.10, pressure altitude from sea level to 
45,000 feet, and free-stream static temperatures to -40° F. 
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SYMBOLS 
The following symbols are used in this report: 
specific heat of iip constant pressure, Btu/(1b)(9F) 
e vapor pressure, 1b/sq ft 
coefficient of evaporation 


ke 
kp coefficient of heat transfer 

L latent heat of vaporization, Btu/1b 
M 


Mach number . 
ma molecular weight of air 
т, molecular weight of water vapor 
р static pressure, Ib/sq ft 
r — factor 


T static temperature, ?R 


t/c airfoil thickness ratio, ratio of airfoil maximum thickness to 
| сһога 


To,e minimum free-stream static temperature corresponding to ice-free 
condition on surface as defined in equation (1а), 


AT- temperature increment resulting from frictional heating in the 
boundary layer, OR А 

Y ratio of specific heats of air, 1-400 

Subscripts: 

О free-stream static conditions 


1 · local conditions at the edge of boundary layer 


8 surface 


= e ee 2 шч. 
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ANALYSIS 


The method of reference 1 offers a convenient means for determining 
the local surface temperature on а wetted'body ав а function of the 
free-stream conditions of temperature, pressure, vapor pressure, and 
velocity, provided the flow field about the body is known. .Ease of cal- 
culation was achieved by a number of simplifying assumptions. The 
restrictions caused by these assumptions will be discussed in relation 
to the general heat balance for an unheated body flying through icing 
conditions. 


Several investigators have made quite complete analyses of the heat 
balance for a body flying through a cloud (references 2 and 5). The 
generalized heat balance for such conditions is well known and has been 
stated quite concisely by Messinger (reference 2). Some simplification 
of the general equation may be obtained by considering the case of an 
insulated body subjected to sufficient frictional heating to maintain 
the surface ice free. In this case, the heat balance is as follows: 


(1) Heat due'to the frictional, or viscous, effects 
plus 
(2) Heat due to the kinetic energy of the water droplets 
plus 
(3) Heat of fusion 
equals 
(4) Heat lost by convection 
plus 
(5) Heat for evaporation of water 
plus 


(6) Heat required to raise temperature of water droplets from stream 
temperature to surface temperature 


The result of Hardy's relation (reference 1) is obtained by equat- 
ing the frictional term (1) to the sum of the convective term (4) and 
the evaporation term (5). For Hardy's relation to be valid, terms (2) 
and (3) must be small compared with term (1), and term (6) must be small 
compared with term (4) or (5); or the sum of terms (2) and (3) must be 
nearly equal to term (6). 


4 RS | МАСА TN 2861 


The heàt of fusion беш (5) ‘offers considerable complication to the 
problem. It їз conceiv&ble that the surface film might consist of a 
mixture of water апа ісе particles without ice actually adhering to the 
surface. The heat, of. fusion would necessarily have to be withdrawn from 
@ portion of the film ‘for such partiéles to occur and, because of the 
lack of knowledge of the percentage of the film that might be frozen 
before visible icing.occurs; it would appear advisable to define an ice- 
free surface as being fully wet: ад, 399 F with no ice particles in the 
filu. This definition.for the icing limit eliminates term (3) from the 
general heat balance and Hardy's relation becomes applicable if term (2) 
is approximately equal Фо екш (6) or both are negligible. Term (2) is 
directly propartional to. the impingement. rate and the square of the 


flight speed; term (6). is directly proportional to the impingement rate . 


and the temperature. difference: between ` the stream and the surface. If 
the impingement rate is small, ‘as for the case of low liquid-water 
contents, then both ‘berms (2) and (6) .&re negligible. For calculations 
of the icing limit; ‘Both terms. are. of the ‘same order of magnitude and 

. because of the nature of the problem the temperature difference between 
the surface and the, stream increases with increasing flight speed. Both 
terms therefore increase with increasing Mach number; large values for 
either term are only encountered ‘at large values of liquid-water content 
and flight speed, Such a combination is highly improbable in actual 
flight since high speeds usually occur at high altitudes, and high alti- 
tudes imply at Jigüld-water contents. 1 


The А of Hardy applies in particular to & surface barely 
wetted with a yery thin film free of ice particles but is probably 
vocis i: correct for не-не капце of practical interest for high-speed 


flight. 


! 


Method of Calculation 
Wetted-surface temperature. - The equation for the temperature of 


a wetted surface in an air stream developed in reference 1 was written 
as 


Sele |. ев e] L | (1) 


where the subscript в refers to conditions at the surface and sub= 
script 1 to conditions at the-edge of the boundary layer. From the 
energy relation for a viscous fluid i 


r(y - 1) „2 
T- eam em [1+ 202) y? 


Substituting this relation in equation (1) апа solving for Tg yields 


POL? 
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Because interest in this problem is fixed оп а definition of flight 
circumstances that provide local surface temperatures of 52° Е, the ` 
terms Ти, І, and eg аге constants. . Hardy observed (reference 1) that 
for the range of temperatures tear 32° К the ratio of the evaporation 
coefficient К. to the heat-transfer ‘coefficient k, was very nearly 
1. A solution of equation (1) can be obtained for any one of the four 
variables (еу, ру, Му, and Ту) provided the ‘other three are known. 


From purely physical reasoning, however, an additional restriction exists. 
This restriction results from the fact that the air in а cloud-is fully 
saturated at the static or free-stream conditions. For each value of 
stream-static temperature, the free-stream vapor pressure is therefore 
assumed a constant and equal to the:saturated vapor pressure. If it is 
assumed that the flow about the body, outside. the boundary layer, is 
accomplished with no change in ‘phase, that is, no condensation or evap~ 
oration, then Dalton's law of partial pressures applies апа 


Pd 


If a reasonable value of tle recovery factor т ів assumed, say 
0.88, then equation (1) may be rewritten in the form shown below 


2 D 0 | 
T. (14 0.176 Mj?) ` x e 
1 И 1. . i 12.75 О 
492 = To - 2776.5 - 
To c ое ~ A/P (Po - eo) 
— pa < 12.75 
Po’ 0 


(1a) 


which is quite convenient for purposes .of calculation because 
Ty/T9 c (1 + 0.176 Mj7^) and pj/Pg can be readily calculated for each 


point on a diamond airfoil. 


Because eg is a function of; То,с, the solution of equation (la) 
must be solved by trial апа error. In an effort to reduce the laborious- 
ness of such calculations, a number of charts similar to that shown in 
figure 1 were made. For each pressure altitude considered, the free- 
stream static temperature corresponding to an ice-free condition To, è 
was plotted as a function of the parameter Т1/Т0,с (1 + 0.176 M) by 
means of equation (la) for constant values of the pressure ratio ру/ро- 
For each value of ру/ро and T4/Tg , (1 + 0.176 М2) а free-stream 
static temperature which corresponds Фо a surface temperature of 520 F 


and therefore represents the minimum free-stream temperature for an ice- 
free surface То, is obtained. The lower limit of free-stream static 


* 
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temperature was considered to be -40° F for the purposes of these cal- 
culations since it has been shown (references 4 and 5) that supercooled 
water droplets are not likely to be found at temperatures below this 
value. 


Flow-fieid calculations. - The flow field about a diamond or double- 
wedge profile can be characterized in general by four regions; a subsonic 
region, a subsonic region with sonic velocity at the shoulder, a mixed 
subsonic-supersonic region with a detached bow wave, and а supersonic 
region with attached shock waves. The regions of most interest are in 
the currently important ranges of transonic and supersonic flight and 
the results presented herein will be limited thereto. 


2704 


А recent analysis is given in references 6 and 7 for the sub- 
sonic region with a sonic line at the shoulder. This solution is for 
a thin, symmetrical wedge at zero angle of attack Pollowed by a straight 
section. The results in this region presented herein are therefore 
limited to the front half of the diamond airfoil. 


The mixed subsonic-supersonic region with detached bow wave has been 
solved in references 8 and 9. This solution applies to a thin, symmet- 
rical diamond airfoil at zero angle of attack for Mach numbers from l to 
the Mach number at which the shock attaches. The flow in this region is 
subsonic over the front half of the airfoil increasing to sonic velocity 
at the shoulder. The flow then expands supersonically around the corner, 
and supersonic flow exists over the back half of the airfoil to the 


trailing edge. 


The flow about a diamond airfoil in purely supersonic flow has been 
known for many years. The results used herein are those of reference 10 
because of their ease in use and applicability. The flow in this region 
is supersonic over the front half of the airfoil with an expansion 
region at the shoulder апа а consequent higher Mach number on the back 
surface of the airfoil, the velocities &nd pressures being constant over 
each surface for а particular value of airfoil thickness and stream Mach 


number. 


The calculation of the chordwise pressure and Mach number distri- 
butions over a symmetrical-diamond airfoil of specific thickness ratio 
t/c at zero angle of attack through the entire transonic-supersonic range 
can therefore be accomplished by combining the results of references 6 


to 10. 


The analyses of references 7 to 9 were made using a speed function f 
and а generalized pressure coefficient Cp: The results of these analyses 
are presented graphically by plots of the chordwise distribution of Ё 
and Cp for constant values of the free-stream speed function EQ. The 
pressure ratio P,/Po is related to the generalized pressure coefficient 


Cp by the equation 
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P 

nonc (аў, 

Po TTE 
2(v +1)5 


The Mach number, either local or free stream, is related to the 
speed function Ё by the equation 


: 2 
1 + |е + (2) 


The temperature ratio T4 /To is obtained from the energy relations 


for & compressible fluid and is given in terms of Mach number by the 
equation 


E 2 
kh в 


"und ede | 


Àn example of & typical calculation of Io e at а particular 
chordwise position is given for the 25-percent chord station and a value 
of the free-stream speed function .Ё of 0.703. Values of Ё of -0.740 
and Cp of 2.890 are obtained from figures 5 and 5, respectively, of 
reference 8. If t/c of 0.06 is assumed, then Mo = 1. 092, М1 = 0.8926, 
ру/ро = 1.276, T,/To = 1.088, and T,/To (1 + 0.176 М1?) = 1.218. At a 
pressure altitude of 15,000 feet, & To, с of 422.5? R is obtained from 
figure 1. 


For the supersonic case with an attached shock (reference 10), the 
values of 51/99 and Му /Mo can be obtained directly from the charts of 


the reference report and, hence, the calculation of the parameters 
required for obtaining То, е from figure 1 is simple. 


E AND DISCUSSION 


For each value of airfoil thiainess ratio and free-stream Mach 


. number assumed, the values of the parameters T,/To (1 + 0.176 My) and 


91/90 can be calculated over the whole range of transonic and supersonic 


speeds. At each pressure altitude;' the’ free-stream static temperature 
To,¢ corresponding to a wetted surface. temperature of 32° F (ice free) 


can then be obtained from charts similar to figure 1; 
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This value of To,c corresponds to the lowest free-stream static 


temperature То which will give ice-free conditions at the particular 
point on the airfoil under consideration for the given values of free- 
stream Mach number and pressure altitude. 


A plot of To, с’ 88 a function of the chordwise distance is shown 


in figure 2 for several free- -stream Mach numbérs for a 4-percent-thick 
symmetrical-diamond QUEE at zero angle of attack and a pressure alti- 
tude of 30,000 feet. 


А study of figure 2 shows that the, general trend of the curves is 
very similar for Mach numbers from O. 7595: to 1.105, that is, from the 
subsónic region to the supersonic region with attached shock. For each 
curve tbe minimum To, c increases slowly from the leading edge to the 


шіасһога. The sudden jump in the free-stream static temperature required 
for icing protection at the midchord results from the expansion around 
the corner and the consequent decrease in local static pressure and 
increase іт local velocity. The increase in velocity and the decrease 

in pressure enhance the evaporation, hence, the increase in the required 
To,c- Water impinging on the front half of the airfoil might therefore 
run back along the surface, spill around the shoulder, and freeze. 
Downstream of the shoulder, the free-stream temperature required for 
protection decreases slowly to the trailing: edge. 


When, the shock wave is attached and,the flow is everywhere super- 
sonic (Mo = 1.2, or greater), the free-stream temperature for protection 
is a constant for the whole front surface and the whole back surface. 

The free-stream temperature required for protection of the back half is, 
however, greater for the reasons previously cited. 


This and similar figures for other altitudes and airfoil thickness 
ratios show that the critical points for these airfoils are immediately 
ahead of the shoulder for the front half of the dirfoil and immediately 
behind the shoulder for the back half of the airfoil. The icing limit 
for these two critical points is shown in figure 3(a) where the minimum 
free-stream temperature corresponding to an ice-free condition is shown 
as a function of the free-stream Mach number for pressure altitudes of 
sea level and 30,000 feet for a 6-percent-thick airfoil. The value of 
То,с decreases rapidly with increasing Mach number; that is, the tem- 


perature region of possible icing hazard. becomes much smaller with 
increasing Mach number. The free-stream static temperatures required 
for protection apparently are higher at the higher pressure altitudes as 
would be expected since, with other conditions equal, the decreased 
pressure enhances the evaporation rate. At a pressure altitude of 
30,000 feet the curve shown for the critical point behind the shoulder 
indicates the possibility of icing at a Mach number approaching 1.4. 
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The effect of altitude is shown clearly in figure 32(b) for free-stream 
Mach numbers of 0.67, 1.0, and 1.5 for the critical point ahead of the 
shoulder on the front surface and for Mach numbers of 1.0 and 1.3 for 
the point immediately behind the shoulder. At a particular value of 
Mach number the curves for both critical points remain nearly parallel 
over the altitude range. For Mach numbers. of 1.0 or less the free- 
stream temperature required for protection increases with increasing 
altitude with a definite decrease in the slope of the -curves in the 
higher altitude range. For the higher Mach numbers the.slope of the 
curves increase over the whole range from sea’ level to 45,000 feet. 


The effect of airfoil thickness on the icing limit for the two 
critical points on the airfoil is shown in figure 3(c) where the free- 
Stream temperature required fer protection is plotted as a function of 
the airfoil thickness ratio.‘ The icing limit at a pressure altitude of 
45,000 feet is shown' for the point ahead of the shoulder for free-stream 
Mach numbers.of 0.875, 1.00, and 1.50 and for the point immediately 
behind the shoulder for Mach numbers of 1.0 and 1.3. For the point 
ahead of the shoulder, the temperature required for protection remains 
constant or decreases with increasing thickness depending on the Mach 
number. : For the point immediately behind the shoulder, the temperature 
required for protection increases with азза thickness. 


An increase in To , of approximately 200 R resuitea when the 


thickness ratio increased from 0.02 to 0-10 for the critical point | 
behind the shoulder at a free-stream Mach number of 1.00. · А decrease in 
To, с of approximately 14.59 R resulted when the thickness ratio 


increased from 0.02 to 0.10 for the critical porn; ahead of the shoulder 
at а free-stream Mach number of 1. 30. 


Experimental verification of the results presented herein has been 
obtained from data taken in a 4- by 10-inch duct tunnel at a free-stream 
Mach number of 1.36 with a 6.82-percent-thick symmetrical-diamond airfoil 
of 5.50-inch chord at zero angle of attack. The icing limit was deter- 
mined by successively raising the free-stream static temperature until 
no ice formed on the surface. The results of these tests аге summarized 
in the following table and compared with results calculated by the 
analytical method presented herein: 


Free-stream | Minimum free-stream static temperature 


static corresponding to ice-free condition, To e 
Pressure Analytical 
Po | 
(in. Hg) Front Back Front Back 
surface surface surface | surface 
9 е 15 ) oda cd „ 
9.41 422 
9.76 TN 
9.91 422 


2.91 
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The excellent agreement obtained in several cases is probably somewhat 
fortuitous since the experimental error was estimated to be 410° F. 
These results, however, show that the calculated icing limit is very 
closely approximated experimentally and, hence, provides considerable 
credence for the calculations presented herein. The run-back icing on 
the rear surface, which was indicated by the analysis, was also observed 
experimentally. Visual observations of the model showed that the surface 
water film flowed smoothly around the shoulder and quickly froze. The 
results presented in the table show that the free-stream static tempera- 
ture at which the run-back freezing occurred experimentally is in good 
agreement with the analytical values. 


The potential hazard to aircraft and missile operation in the tran- 
sonic region may be seen by applying the results previously given for 
the critical points at the shoulder of a diamond airfoil to operation 
under atmospheric conditions likely to be encountered: If a particular 
relation of pressure and temperature with altitude in the atmosphere is 
assumed, it is possible to calculate the flight Mach number necessary 
for icing protection for a given airfoil thickness as a function of 
altitude. The Mach number necessary for protection for the critical 
point ahead of the shoulder for a 6-percent-thick diamond airfoil is 
shown in figure 4 for NACA standard atmosphere (reference 11). For these 
conditions the region of icing hazard exists from approximately 7500 to 
27,300 feet. If a more severe temperature relation with altitude is used 
such as might be encountered during winter, that is, a 52° F sea-level 
surface temperature: and a lapse rate of 3° F per thousand feet, then the 
curves shown result for the critical points at the shoulder. For these 
conditions the region of hazard exists from sea level to 24,000 feet. 

A Му of 1.29 is needed at the peak altitude in order to provide 


protection. 


Although the results presented in this report apply only to the 
symmetrical-diamond airfoil at zero angle of attack, they may be used 
as an indication of the region of icing hazard and hence serve as a 
guide for determining where de-icing systems are needed or what flight 
plan may be used to prevent icing. 


SUMMARY OF RESULTS 


The results of the analysis of icing limit presented for the 
diamond airfoil at zero angle of attack in the transonic-supersonic 
speed range are ав follows: 


1. The critical point for the whole airfoil in terms of the icing 
limit was located immediately behind the shoulder. The critical point 
on the front surface was located immediately ahead of the shoulder. 
The results indicate that icing may occur at the critical point behind 
the shoulder for Mach numbers as high as 1.4. 
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2. At а constant free-stream Mach number апа airfoil thickness, the 
minimum free-stream static temperature corresponding to an ice-free 
condition increases with аг: altitude, 


5. For the point ТР Та ahead of the shoulder the minimum 
free-stream temperature corresponding to an ice-free condition remains 
the same or decreases with increasing thickness ratio at & constant 
value of Mach number and altitude. For the point imtediately down- 
Stream of the shoulder the free-stream temperature corresponding to 
an ice-free condition increases with increasing thickness at a constant 
Mach number and altitude. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 24, 1952 
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Airfoil thickness, 4 percent; pressure altitude, 50,000 feet, 


for several Mach numbers, 


Figure 2 
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(в) Effect of Mach PIEP Pressure altitudes, sea level and 30,000; 
airfoil thickness, 8 percent. 
Figure 3, - Variation of minimum free-stream static temperature corresponding to ice-free condi- 
tion for two critical points on airfoil. 
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(b) Effect of pressure altitude for several Mach numbers, Airfoil thickness, 6 percent. 


Figure 3. - Continued. Variation of minimum free-atream static temperature corresponding to ice-free condition 
for two critical points on airfoil. 
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(c) Effect of airfoil thickness for several Mach numbers. Pressure &ltitude, 45,000 feet. 


Figure 5. - Concluded. Variation of minimum free-stream static temperature corresponding to ice-free 
condition nd two critical points on airfoil. 
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Figure 4. - Variation of flight Mach number required for icing protection as 
function of altitude for several assumed atmospheric temperature-pressure 
relations.  Airfoil thickness, 6 percent. 
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